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THE  COMPONENTS  OF  POWER  APPEARING  IN  THE 
HARMONIC  ANALYSIS  OF  A STATIONARY  PROCESS 


by 


M.  M.  Siddiqui 


1.  INTRODUCTION 


Let  Y(t)  be  a real  stationary  process  with  mean  zero  and  auto- 
covariance  function  C(t)  = EY(s)  Y(s+t).  Here  Ex  denotes  the  expected 
value  of  the  random  variable  x and  we  assume  C(t)  to  be  finite  for  all 
values  of  t.  Suppose  we  have  a continuous  record  of  Y(t)  over  the 
time  interval  (0,T)  measured  in  seconds  and  we  take  2n+l  equally 
spaced  readings  YQ,  Yj,  ...»  Y£n  at  times  t=0,  6,  26,  . . . , 2n6  = T 
seconds.  It  is  customary  to  represent  the  readings  by  the  Fourier 
aeries 


+ {3j  sin 


» • • • » 


to  estimate  the  power  spectrum  of  the  process  Y(t). 


Assuming  that  the  true  representation  of  Y(t)  is  given  by 


oo 


(1.2)  Y(t)  = S 

j=l 


1 


it  will  be  shown  that  the  power  ascribed  to  the  frequency  j/T  cycles 
per  second  (cps)  in  the  analysis  of  (1, 1)  consists  of  three  components: 

(1)  the  true  power  in  the  frequency  j/T  cps; 

(2)  the  powers  in  the  frequencies  (2kn  + j)/T  cps,  k=  1,  2, . . . ; 

(3)  a part  of  the  powers  in  the  frequencies  l/(mT)  cps, 
m 2,  3,  •••• 

The  frequencies  mentioned  in  (2)  are  called  the  Aliases*  to 
the  frequency  j/T  cps  and  the  transfer  of  their  powers  to  the  frequency 
j/T  cps  is  well  known.  The  study  of  the  effect  of  the  presence  of  fre- 
quencies lower  than  l/T  cps  in  the  process  on  the  harmonic  analysis, 
however,  seems  to  have  been  somewhat  neglected.  Spetner  [1]  con- 
siders the  error  in  power  spectra  due  to  finite  sample  from  a con- 
tinuous Gaussian  process,  and  to  the  author* s knowledge  no  other 
references  are  found  in  the  literature.  The  following  results  do  not 
suppose  the  process  to  be  Gaussian  and  the  only  restrictions  placed  on 
the  process  are  given  by  equations  (2.  2)  and  (2.  3). 

2.  LOW  FREQUENCIES  IN  THE  PROCESS 

The  record  of  Y(t)  over  the  interval  of  time  (0,T)  can  be 
exactly  represented  by 

(2.  1)  Y(t)  = ac  + 2 (a.  cos  + pj  sin 

j=l  V 


2 


However,  let  the  true  representation  of  the  process  Y(t)  be 


(2.2) 


Y(t)  = 


oo 
+ S 
m=2 


cos  + P*(j)  sin 

(a*(i)  cos  gr  + P*  <^>  ^ 


2irt 

mT 


where  a*  and  (3*  are  random  variables  satisfying  the  following  con- 
ditions: 


(2.  3)  Ecl*(s)  = E{3*(s)  = 0 for  all  index  values  of  s , 

Ea*2(s)  = Ep*2(s)  = <r*2(s)  , 

Ea*(s)a*(r)  = Ep*(s)|3*(r)  =0  if  s ^ r , 

Ea*(s)p*(r)  = 0 for  all  s and  r , 

E o-*2(s)  = (Ty2  < oo  . 
s 

The  power  in  the  frequency  s/T  cps  is  then  defined  as  cr*2(s)  . 

The  conditions  (2.  3)  are  satisfied,  for  example,  if  the  typical 

term  of  (2.  2)  is  represented  by  yg  cos  ^ 2]lst,  + c8  ^ and  we  assume 

2 2 

yg  and  cg  to  be  independent  variates  with  Eyg  = 0 , Ey  s = 2o-*  (s) 
and  € to  be  distributed  uniformly  over  the  interval  (0,2tt). 

Now,  for  0 < t < T 
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CD 


(2-4) 


2irt  rn.2irm.2ir  _ ,.Z  2 1X-1  2irjt 

cos  — — • = — — sm Bin  — S (i  m -1)  cos — 

mT  2ir  m ir  m . 1 T 


7 oo 

, mc  n 2*n\  — .,.2  2 , 2irjt 

+ (1  - cos  — ) S i(j  m -1)  sin  ■ * 

ir  m T 


CD 

2irt  m /,  2tt.  m n 2ir.  — «-  „ 

sin  — — = — (1  - cos  — ) - — (1  - cos  — ) S (j  m -1)  cos 
ml  £>tt  m ir  m j_  ^ 1 


2 2 ,rl  .o2^ 

I m -1 » COS  ■ — * — 


OO 


^ sin  £ S j(j2m2-l)_1  sin  2^1 

ii  m j=1  T 


Since  (2.1)  and  (2.2)  are  identical  over  the  interval  (0,T), 


we  obtain 
(2. 5)  a_ 


= s ro*(— )|2- sin  — + p*(— )2L  (1  - cos  — )1  . 
J m 2ir  m m 2ir  m l 

m=2u  J 


and  for  j > 0 


oo 


(2.6)  a.5  = a*(j)  - S 

m=2 


a*(— )(j^m^-l)“*  — sin  — + p*(—  )(j^m^-l)~* 
m ir  m m 


m 2ir  , I 

— (1  - cos  — ) 

* m J 


Pi  = P*0)  + S f a*(-L)  j(j2m2-lf 1 m2(l  . cos  i!)  - p*(-L)j(j2m2-lf 1 

._->L  m IT  m 


m=2 


m 


“-sin2! 


7T 


m 


]• 


We,  therefore,  have 


(2.7)  Ea.j  = EPj  = 0 for  j = 0,1, ..., 
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CD 


(2.8)  Ea_  = 2 <r*  (— ) 


2,  1 , m2  2 ir 


m=2 


m n2 


sin 


m 


E.f.,.2U)  + 4 ” "2  •‘°2  g 

J m=2  m ir2(j2m2-l)2 


5 , °°  ? 1 m^j2  sin2  — 

Epj  = <r*2(j)  + 4 2 <T*2(-i-)  ™ 

J m=2  m ir2(j2m2-l)2 


00  2 1 m2  sin2  wr 

Ea0aj  = -2  2 cr*' 2(^-> 

J m=2  m ir2(j2m2-l) 


00 


r2  -2  TT 

m sin  — 
m 


Ea^as  =42  cr*2(J-)  - , 

m=2  m ir2(j2m2-l)(i2m2-l) 


OO  ij  sin2  JL 

Ep.p.  = 4 S o-^(— ) no 

1 J -o  m 2/. 2 2 1W,2  2 

m=2  tt  (j  m -l)(i  m -1) 


for  j > 0 , 


for  j > 0 , 


for  j > 0 , 


for  i^j,  i ,j / 0, 


for  i / j , 


Eojpj  = 0 for  all  i and  j . 


W riting 


o2  = Eoq2,  o-j2  = | E(  cij2  + Pj2  ) for  j > 0 


we  find 


(2.  9)  <r  2 = 2 f- sin  — Y <r*2  (JL)  , 

m=2  ^ * m ' m 


, 2 00  j2m2+l 

of  = <r*  U)  + 22 
•/ 


2 2 ix2 


m=2  (j  m -1) 


(—  sin  — ^ <r*2(— ),  for  j >0, 
\ it  m / m 
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Performing  the  calculations  so  as  to  account  for  at  least  90% 
of  the  power  in  each  frequency,  we  have 

(2.10)  0-o2  = 0.405  <r*2(I)  + 0.684  <r*2(I)  + 0.  811<r*2(I)  + 0.  874<r*2(|-) 

+ 0.  912  <r*2(i)  + 0.933  <r*2d)  + 0.950  o-*2(I)  + e S o-*2(— ) 

6 7 8 m=9  m 

where  0.  96  < 0 < 1 , 

O-  2 = <r*2(  1)  + 0.450  <r*2(— ) + 0.  214  (r*2(!)  + 0.123  <r*2(I) 

i i 2 3 4 

+ 0.  079  <r*2(h  + . . . , 

D 

o>2  = <r*2(2)  + 0.  061  + 0.  041  o-*2(i)  + . . . , 

a Z 3 

o-32  = <7*2(3)  + 0.  024  o-*2(j)  + 0.  018  o-*2(-|)  + . . . . 

Oj2  S <r*2(j)  for  j > 4 . 

We  observe  that  more  than  95%  of  the  power  in  the  frequencies  lower 

than  l/(7T)  cps  is  transferred  to  the  zero  frequency.  Since  the  mean 

of  the  sample  is  usually  subtracted  from  each  observation  before  fur- 

2 

ther  analysis,  the  power  o*Q  will  be  subtracted  from  the  spectrum. 

Furthermore,  more  than  90%  of  the  power  in  frequencies  ranging 

from  — cps  to  zero  cps  is  transferred  to  the  zero  and  the  first  two 
2T 

lowest  frequencies  analysed. 
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The  scope  of  the  scientific  program  of  the  National  Bureau  of  Standards  at  laboratory 
centers  in  Washington,  D.  C.,  and  Boulder,  Colorado,  is  given  in  the  following  outline: 

Washington,  D.C. 

Electricity  and  Electronics.  Resistance  and  Reactance.  Electron  Tubes.  Electrical  Instruments. 
Magnetic  Measurements.  Dielectrics.  Engineering  Electronics.  Electronic  Instrumentation. 
Electrochemistry. 

Optics  and  Metrology.  Photometry  and  Colorimetry.  Optical  Instruments.  Photographic 
Technology.  Length.  Engineering  Metrology. 

Heat  and  Power.  Temperature  Physics.  Thermodynamics.  Cryogenic  Physics.  Rheology  and 
Lubrication.  Engine  Fuels. 

Atomic  and  Radiation  Physics.  Spectroscopy.  Radiometry.  Mass  Spectrometry.  Solid  State 
Physics.  Electron  Physics.  Atomic  Physics.  Nuclear  Physics.  Radioactivity.  X-rays.  Betatron. 
Nucleonic  Instrumentation.  Radiological  Equipment.  AEC  Radiation  Instruments. 

Chemistry.  Organic  Coatings.  Surface  Chemistry.  Organic  Chemistry.  Analytical  Chemistry. 
Inorganic  Chemistry.  Electrodeposition.  Gas  Chemistry.  Physical  Chemistry.  Thermo- 
chemistry. Spectrochemistry.  Pure  Substances. 

Mechanics.  Sound.  Mechanical  Instruments.  Fluid  Mechanics.  Engineering  Mechanics. 
Mass  and  Scale.  Capacity,  Density,  and  Fluid  Meters.  Combustion  Controls. 

Organic  and  Fibrous  Materials.  Rubber.  Textiles.  Paper.  Leather.  Testing  and  Specifica- 
tions. Polymer  Structure.  Organic  Plastics.  Dental  Research. 

Metallurgy.  Thermal  Metallurgy.  Chemical  Metallurgy.  Mechanical  Metallurgy.  Corrosion. 
Metal  Physics. 

Mineral  Products.  Engineering  Ceramics.  Glass.  Refractories.  Enameled  Metals.  Concreting 
Materials.  Constitution  and  Microstructure. 

Building  Technology.  Structural  Engineering.  Fire  Protection.  Air  Conditioning.  Heating,  and 
Refrigeration.  Floor,  Roof,  and  Wall  Coverings.  Codes  and  Specifications.  Heat  Transfer. 
Applied  Mathematics.  Numerical  Analysis.  Computation.  Statistical  Engineering.  Mathemat- 
ical Physics. 

Data  Processing  Systems.  SEAC  Engineering  Group.  Components  and  Techniques.  Digital 
Circuitry.  Digital  Systems.  Analogue  Systems.  Application  Engineering. 

• Office  of  Basic  Instrumentation  • Office  of  Weights  and  Measures 

Boulder,  Colorado 
BOULDER  LABORATORIES 
F.  W.  Brown,  Director 

Cryogenic  Engineering.  Cryogenic  Equipment.  Cryogenic  Processes.  Properties  of  Materials. 
Gas  Liquefaction. 

Radio  Propagation  Physics.  Upper  Atmosphere  Research.  Ionospheric  Research.  Regular 
Propagation  Services.  Sun-Earth  Relationships. 

Radio  Propagation  Engineering.  Data  Reduction  Instrumentation.  Modulation  Systems.  Navi- 
gation Systems.  Radio  Noise.  Tropospheric  Measurements.  Tropospheric  Analysis.  Radio 
Systems  Application  Engineering. 

Radio  Standards.  High  Frequency  Electrical  Standards.  Radio  Broadcast  Service.  High  Fre- 
quency Impedance  Standards.  Calibration  Center.  Microwave  Physics.  Microwave  Circuit 
Standards. 
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